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ABSTRACT 


A  satellite  has  been  observed  through  an  optical  telescope 
with  a  CCD  camera,  the  electrical  signal  stored,  and  later 
analyzed  with  the  aid  of  a  digital  computer.  The  unresolved 
satellite  image  is  contained  in  3  to  7  picture  elements  in  fifty 
contiguous  data  frames.  This  allows  a  better  understanding  of 
the  physical  characteristics  of  the  satellite  image  as  viewed 
through  the  atmosphere  by  an  electro-optical  detector. 

This  paper  describes  the  data  processing  programs  used  to 
analyze  the  data,  the  analysis  results,  and  implications  for 
signal  processing  design. 
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I .  INTRODUCTION 

A  CCD  camera  has  the  potential  of  improving  the  performance 
of  satellite  detection  systems  by  means  of  higher  quantum 
efficiency,  smaller  pixel  size,  and  more  pixels  on  a  focal  plane 
than  EBSICON  cameras.  An  efficient  satellite  detection  processor 
should  be  designed  to  exploit  the  signal  characteristics  of  the 
camera  it  uses  as  an  electro-optical  detector.  Selected  portions 
of  a  digital  satellite  detection  processor  (IMC  Signal  Processor) 
and  CCD  camera  are  being  developed  for  the  GEODSS  advanced  tech¬ 
nology  program.  The  output  of  a  CCD  camera,  when  observing  a 
satellite,  has  been  digitized,  and  then  analyzed  with  the  aid  of 
a  computer.  This  allowed  a  detailed  study  of  the  electrical  signal 
that  the  processor  would  have  to  use  prior  to  designing  the 
processor. 

Satellite  ATS-6  was  observed  with  the  31"  main  GEODSS  ETS 
telescope  at  the  White  Sands  Missile  Range.  The  CCD  camera  was 
at  the  prime  focus  and  was  operated  in  the  stare  mode  with  an  0.1 
second  integration  time.  The  ETS  photometer  measured  the 
satellite's  brightness  as  10.5  Mv.  The  camera  video  was 
digitized  to  12  bit  resolution,  sent  to  the  MODCOMP  computer 
at  the  ETS  and  stored  on  magnetic  tape.  This  tape  was  brought 
back  to  Lexington  and  placed  on  the  MODCOMP  computer  here.  The 
tape  was  read  by  the  MODCOMP  and  the  data  transmitted  to  the 
HP  9845  graphics  computer  and  stored  there.  The  HP  computer 
was  used  for  all  the  data  analysis. 


The  satellite  observation  is  contained  in  fifty  contiguous 
frames.  Each  frame  contains  an  array  of  32  x  32  picture  elements 
(pixels).  The  integration  time  for  each  frame  was  0.1  second. 
Each  pixel  has  the  equivalent  field-of-view  of  2" 7  x  2’.’7. 
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II.  DATA  PROCESSING 

The  frames  of  data  were  read  from  the  digital  tape  by  the 

MODCOMP  computer  and  transmitted,  on  an  RS-232C  interface1,  to 

2 

the  HP  computer.  Using  this  interface  the  HP  computer  was 
operated  to  look  like  a  standard  terminal  to  the  MODCOMP.  The 
major  task  of  this  interface  was  to  write  an  HP  terminal 
emulator  program.  This  program  has  to  read  the  data  from  the 
interface,  store  the  data,  and  respond  to  the  MODCOMP  with 
standard  terminal  handshake  commands.  The  program,  "MOD-HP", 
listed  in  Appendix  A  did  this  at  9600  baud. 

To  operate  at  this  high  a  rate,  each  data  word  was  read 
as  a  character  string.  The  character  strings  were  stored  on 
HP  tape  cassettes  as  the  data  word  was  being  read.  The  1024 
character  strings  from  each  frame  were  separated  by  markers. 

Once  all  fifty  frames  of  the  character  strings  were  stored  on 
tape,  the  HP  computer  was  disconnected  from  the  MODCOMP,  and 
was  operated  independently  from  then  on. 

Two  tape  cassettes  were  required  to  store  all  the  character 
strings.  All  the  character  strings  were  then  read  off  the  tapes 
and  stored  on  one  floppy  disc  using  the  "COPY"  program  listed 
in  Appendix  A.  The  character  strings  cannot  be  used  in 
programs,  but  must  be  converted  to  the  numbers  they  represent. 
The  program  "READ",  Appendix  A,  read  the  character  strings  off 
the  disc,  converted  the  character  strings  to  numbers  and  stored 

IT!  Hewlett  Packard,  "Serial  I/O  Interface-98036A, "  1977. 

2.  Hewlett  Packard,  "I/O  ROM  Programming  System  45B,"  1979. 


these  on  tape.  The  "COPY"  program  was  again  used  to  transfer 
data  from  tape  to  disc. 

To  insure  that  the  data  stored  on  the  disc  was  accurate, 
a  print-out  of  the  data  stored  on  the  MODCOMP  tapes  was  compared 
to  the  data  stored  on  the  HP  disc.  The  data  was  identical  in 
frame- to- frame  and  entry-to-entry  comparisons.  The  first  data 
entry  in  each  frame  of  both  the  MODCOMP  and  HP  data  was  the 
same,  but  was  not  a  valid  data  value.  To  avoid  problems  in 
data  processing  the  value  of  the  first  data  entry,  in  each  frame, 
was  set  to  the  mean  of/  the  data  of  each  frame  and  is  then  called 
a  valid  data  frame. 

Now  the  data  is  ready  to  be  analyzed.  The  first  data 
analysis  program  is  "MAXAMP"  listed  in  Appendix  B.  "MAXAMP”  finds 
the  maximum  amplitude  in  each  frame  and  prints  the  frame  number, 
amplitude,  and  array  position.  A  sample  of  "MAXAMP"  output  is 
shown  in  Figure  1.  The  frame  number,  e.g.  T11V01:F  stands  for 
CCD- DIGITIZER  test  #11  (Til) ,  valid  data  frame  (V) ,  frame  number 
(01) ,  and  stored  on  floppy  disc  (F) .  The  pixel  array  position 
is  numbered  from  1  to  1024  in  each  data  frame. 

The  next  data  analysis  program  listed  in  Appendix  B  is 
"FREDIS"  which  plots  the  frequency  distribution  of  amplitudes 
for  each  frame.  Figure  2  contains  plots  from  two  of  the  fifty 
frames.  The  amplitude  scale  adjusts  for  each  frame  depending  on 
the  maximum  amplitude.  The  occurrence  scale  also  adjusts 
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depending  on  the  number  of  pixels,  excluding  the  zero  amplitude, 
for  which  there  are  from  350-400  pixels  in  each  frame. 

The  program  "PEAKS"  examines  a  frame  of  data  and  prints  out 
the  location  and  amplitude  of  any  pixel  whose  amplitude  is  above 
a  selected  level.  Figure  3  is  a  sample  of  "PEAKS"  output.  A 
program  listing  of  "PEAKS"  is  in  Appendix  B.  The  frame  number 
is  given  and  the  selected  level  (e.g.  19).  The  row,  column, 
array  position  (as  in  "MAXAMP"  output) ,  and  corresponding 
amplitude  are  listed. 

Program  "DATAST",  Appendix  B,  prints  out  the  amplitude 
values  of  any  frame  by  row  and  column.  An  asterisk  is  placed 
next  to  amplitudes  above  a  selected  amplitude.  The  output  is 
displayed  on  the  CRT  or  printed  on  the  line  printer.  A  sample  of 
output  from  frame  #11  is  shown  in  Figure  4.  The  top  line  of 
each  pair  is  array  elements  1  to  16  and  the  second  line  is 
array  elements  17  to  32.  In  this  frame  the  satellite  pixels 
have  amplitudes  of  182,  112,  32,  and  40.  An  isolated  pixel  of 
noise  with  amplitude  32  is  one  array  line  above  this. 

The  "FRAME"  program  (Appendix  B)  allows  a  13  x  11  pixel 
array  around  the  satellite  image  of  any  frame  to  be  displayed. 
After  the  array  of  pixels  is  displayed,  the  amplitudes  are 
plotted  first  in  column  and  then  in  row.  Figure  5  is  a  print¬ 
out  of  the  array  of  pixels  around  the  satellite  image  from 
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frame  #12.  Figure  6  is  a  plot  of  the  amplitude  in  the  sixth 
column  of  Figure  5.  The  dashed  line  is  a  level  at  an  amplitude 
of  19.  Figure  7  is  a  plot  of  amplitudes  in  the  fifth  row  of 
Figure  5.  The  dashed  line  is  also  at  an  amplitude  of  19. 

Appendix  C  contains  the  pixels  (5  x  6)  surrounding  the 
satellite  image  for  frames  1  through  12.  These  were  obtained 
using  a  modified  version  of  program  "FRAME" . 


III.  ANALYSIS  OF  DATA 

Analysis  of  the  output  of  the  "MAXAMP"  program  indicates 
that  the  highest  amplitude  of  the  satellite  image  moves  between 
two  adjacent  pixels.  Over  the  fifty  frames,  the  highest 
amplitude  switches  15  times  between  the  two  pixels,  spends  an 
equal  amount  of  time  in  each  and  there  is  no  periodicity  to  the 
switches.  The  range  of  highest  amplitudes  is  127  to  224  with 
a  mean  of  165.  The  mean  of  highest  amplitude  in  one  pixel  is 
162,  while  the  other  pixel's  highest  amplitude  mean  is  169. 

Examining  the  output  of  "FREDIS"  shows  that  in  all  frames 
most  of  the  pixels  have  amplitudes  from  0  to  17.  The  mean 
of  the  amplitudes  from  each  of  fifty  frames  ranges  from  4.4  to 
4.9.  This  is  strongly  due  to  the  large  number  of  pixels  (350  to 
400)  in  each  frame  with  amplitudes  of  zero.  The  distributions 
of  all  fifty  frames  are  similar.  All  frames  have  less  than  20 
pixels  above  an  amplitude  of  19. 

The  output  of  "PEAKS"  shows  that  in  every  frame  the  pixels 
with  amplitudes  greater  than  28  are  either  those  of  the  satellite 
image  or  of  isolated  pixels.  There  are  532  pixels  that  have 
amplitudes  greater  than  19  in  the  fifty  frames.  Most  of  these 
pixels  are  either  those  of  the  satellite  image  or  isolated 
pixels.  However,  12  of  the  pixels  in  10  of  the  fifty  frames 
have  an  adjacent  pixel  greater  than  19.  In  half  of  these  cases 
the  adjacent  pixels  are  vertically  related  and  in  half  horizontally. 
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The  pixels  with  amplitudes  greater  than  39  are  all  those  of 
the  satellite  image. 

The  output  of  "DATAST"  shows  that  the  noise  is  uncorrelated 
from  pixel  to  pixel  and  from  line  to  line.  Eight-four  percent 
of  the  pixels  with  an  amplitude  greater  than  8  have  no  adjacent 
pixel  with  amplitude  greater  than  8.  Ninety-three  percent  of 
the  pixels  with  amplitude  greater  than  12  have  no  adjacent  pixel 
with  amplitude  greater  than  12. 

The  output  of  "FRAME"  gives  information  primarily  about  the 
pixels  associated  with  the  satellite  image.  Table  1  is  a  summary 
of  the  distribution  of  satellite  image  pixels  by  frames.  The 
distribution  of  energy  from  the  satellite  image  is  very  dynamic. 
The  distribution  changes  significantly  in  periods  of  0.1  sec  to 
0.3  sec.  In  61.5%  of  the  frames  the  energy  distribution  doesn't 
fit  a  circular  Gaussian  function. 


IV.  SIGNAL  PROCESSING  CONSIDERATIONS 

After  reviewing  the  analysis  results,  it  is  apparent  that 
the  satellite  image  is  not  stable  over  time  periods  as  short 
as  0.2  sec.  The  instability  of  the  image  could  result  from  three 
primary  sources:  satellite  motion,  telescope  jitter,  and  atmospheric 
seeing.  During  the  time  of  data  collection,  the  ATS-6  satellite 
was  a  three-axis  stabilized  geosynchronous  satellite.  Its 
reflected  light  vs.  time  had  been  measured  by  several  observers 
at  the  ETS  and  was  found  to  be  steady  throughout  the  night,  as 
would  be  expected  from  this  type  of  satellite.  Therefore  satellite 
motion  didn't  account  for  the  image  instability  observed  in  the 
data.  The  fifty  frames  of  data  were  carefully  scrutinized  to 
detect  any  indication  of  telescope  jitter.  Telescope  jitter 
would  be  detected  in  the  data,  if  the  pixels  which  contain  the 
satellite  image  shifted  together  in  the  same  direction.  There  is 
no  observable  correlation  between  the  maximum  amplitude  pixel  shift 
direction  and  the  other  satellite  image  pixels'  shift  direction. 
Calculating  the  centroid  of  the  satellite  pixels  shows  that  the 
centroid  movement  is  in  random  directions  rather  than  in  one  or 
two  orthogonal  directions  that  would  characterize  telescope  jitter. 
Atmospheric  seeing  is  therefore  the  leading  candidate  as  the 
source  of  the  image  instability.  Although  few  measurements  of 
atmospheric  seeing  with  time  constants  as  short  as  0. 1-0.2  sec, 
and  with  two-dimensional  spatial  resolution  as  fine  as  2'.’7  have 
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3 

been  reported,  one  reference  does  report  the  characteristic 
image  fluctuation  seen  in  this  data  as  atmospheric  seeing.  Many 
other  reports  on  atmospheric  seeing  describe  image  fluctuation 
similar  to  that  observed  in  this  data.  Atmospheric  seeing  will 
affect  stellar  images  in  the  same  fashion.  A  signal  processing 
system  that  removes  stars  by  comparing  a  star  field  at  one  time 
with  the  same  star  field  at  a  time  0.2  sec  later  or  longer,  must 
incorporate  a  technique  to  reduce  false  alarms  due  to  the  star 
image  instability.  These  false  alarms  are  located  adjacent  to 
star  images,  and  a  guard  band  around  each  star  that  inhibits 
detections  would  remove  these  false  alarms. 

Although  the  resolution  of  the  telescope  is  1'.'0  or  smaller 
and  the  pixel  resolution  is  2'.'74,  the  satellite  image  is  in 
3  to  7  pixels.  Most  of  the  energy  is  in  2  to  3  pixels,  however 
in  74%  of  the  frames  the  energy  is  in  4  to  5  pixels.  In  more 
than  half  the  frames  the  distribution  of  satellite  energy 
doesn't  fit  into  a  circular  function.  There  is  little  correla¬ 
tion  from  frame- to- frame  of  the  spatial  distribution  of  satellite 
energy.  This  requires  that  predictions  of  satellite  SNR  and 
system  detection  performance  be  modified  to  incorporate  these 
physical  realities. 


3.  MIKESELL,  A.H.  et  al,  "The  Scintillation  of  Starlight",  JOSA, 
41,  689-695,  1951. 
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The  star  and  satellite  images  are  in  several  adjacent  pixels 
(at  a  2'.'74  pixel  size),  while  the  noise  is  only  one  pixel  long. 
This  is  an  excellent  discriminant  between  the  satellite  and  high 
amplitude  noise  that  competes  with  the  satellite  signal.  If 
a  threshold  level  of  28  was  used  with  the  pixels  in  the  fifty 
frames,  the  threshold  crossings  would  be  due  to  the  satellite 
(adjacent  threshold  crossings)  or  to  the  noise  (isolated  thres¬ 
hold  crossings).  The  isolated  threshold  crossings  could  then  be 
removed  and  only  the  satellite  threshold  crossings  would  remain. 
As  the  threshold  is  lowered,  noise  threshold  crossings  will  be 
more  frequent  and  false  alarms  due  to  noise  will  start  to  re¬ 
appear.  At  a  threshold  of  28,  100%  of  the  noise  false  alarms 
are  removed;  at  a  threshold  of  12,  93%;  and  at  a  threshold  of 
8,  84%.  The  remaining  false  alarms  due  to  noise  can  be  removed 
by  additional  processing  (streak,  slope,  velocity,  etc.). 
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V. 


CONCLUSIONS 


The  analysis  of  a  satellite  image  as  seen  by  a  CCD  camera 
mounted  on  the  31"  prime  focus  at  the  ETS  has  been  accomplished. 
The  complexity  of  the  image  as  a  function  of  spatial  distribution 
and  time  can  be  clearly  seen.  The  seeing  effects  were  worse  than 
what  had  previously  been  expected  at  the  ETS.  The  utility  of 
having  a  digitized,  stored  signal  from  a  sensor  before  designing 
a  signal  processor  that  uses  that  particular  sensor,  is  obvious. 

The  pixel  size  of  a  camera  is  important  to  the  detection 
system;  as  too  small  a  pixel  will  cause  the  satellite  energy  to 
be  spread  over  many  pixels  and  thus  reduce  the  SNR,  or  if  the 
pixel  is  too  big  the  night  sky  will  contribute  additional  noise 
and  also  lower  the  SNR.  The  image  size  depends  on  the  atmospheric 
seeing,  therefore,  the  selection  of  pixel  size  depends  on  the 
seeing  expected  at  a  site.  The  ideal  way  to  measure  seeing  at 
the  ETS  would  be  to  use  an  automated  version  of  what  was  done 
in  this  report.  The  microprocessor  based  PVDS  system  would  be 
ideal  for  this  purpose  and  could  provide  a  short  exposure 
measure  of  seeing.  Using  a  CCD  camera  mounted  on  the  cassegrain 
position  of  the  31"  ETS  telescope  would  provide  even  more  detail 
of  the  image  intensity  distribution.  Binary  stars  could  be 
viewed  to  eliminate  some  observing  uncertainties  and  data  could 
be  collected  on  a  regular  basis. 
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1>I 

DIM  D*' 500  '[  3621 

!  "MOD-HP" 

20 

lREhTE  "TtlFOl", 1, 1O00O 

30 

RSSIGN  #1  TO  "T11F01" 

40 

RESET  10 

SO 

WRITE  10  10.5;  1 

15  0 

WRIT  WRITE  10,4564 

70 

WRIT  WRITE  lO. 4; 203 

30 

WRIT  WRITE  10,4; 32+7 

90' 

READ  10  1 0 , 4 ; R  4 

100 

WRITE  10  10, 5; O 

1  1  u 

WRIT  RERI'  10, 4;  R 

1  RERDY 

1  20 

WRITE  10  lO. 7; O 

i TRIGGER 

1  30 

i 

1  4  0 

i  SET  UP  THE  9S0 

36R  TO  INTERRUPT  ON  INPUT 

DRTR 

ISO 

! 

1  6  0 

WRITE  10  10,5; 128+4 

1  70 

P  = I 0  F  L  R  G  >  10  • 

1  30 

8=  1  OSTATUSt  10  > 

1  90 

PRINTER  IS  16 

2'  0  O' 

PRINT  "10  FLRG  1  ",F,"I 

0  STRTUS  1  “  ,  S 

2  1  0 

RERI'  10  10,5;  R5 

2  20 

READ  1 0  1  0 ,  6 ;  P  6 

2  3  O' 

STRTUS  10; V 

24  0 

PRINT  "  R:4  O  "  ,  R4 ,  "  R5 

146  “,R5,"P6  248  ",R6, 

"  V 

4  02  ",  V 

2  S  O' 

1  =0 

! INDEX  FOR 

WORDS 

COMING  OVER 

2  6  O 

BEEF 

270 

BEEF 

230 

TOPEN  10,2  GOSLiB  Get  input 

2  9  0 

Loop: 

1  WRIT  FOR 

I'RT  R 

FROM  MODCOMP 

3  0  0 

ON  KBD  1  GOTO  310 

3  1  0 

K  f  =  K  G  D  t 

320 

E  x  =  P  0  $ «:  K  $ ,  c  H  R  *  '  2  5  5  >  «,  C  H  R  *  <  2 

1  >  > 

3  3  O' 

IF  E  •.  0  THEN  440 

3  4  0 

GOTO  Loop 

350 

Get i n  p  u  t : 

!  GET  INPUT 

FROM 

INTERFACE 

360 

1  =  1  +  1 

370 

IF  I  500  THEN  Loop 

380 

UK  I  >  =TBUF ♦ 

3  9  0 

I  F  POS  <.  B t  ■  I  >  ,  CHR*  <.  5  >  >  THEN 

OUTPUT  10; CHR# C 6 > 

4  0  0 

PRINTER  IS  16 

410 

PRINT  I 

420 

WRIT  90O 

4  30 

RETURN 

440 

FOR  1=1  TO  500 

450 

PRINT  tl  ;  I  ,  Ii*<  I  > 

460 

NEXT  I 

470 

END 

"MOD-HP" 
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FOR  k-J  =  1  0  TO  50  !  "feEAD" 

D  I  H  D*c2O0>  C  320]  ,  PC 321  > ,  0<  200 >  .  F  <  1  024  :> 
x#=val*‘ w> 

U I  *  =  "  T 1  IF  "8.X*?  “:F8" 

W2#  =  "  T  1  1  D  "  &X#&  "  :  T15" 

ASSIGN  #1  TO  Ml* 

ON  END  #1  GOTO  680 
St  =0 
D  =  0 

FOR  1^1  TO  20O  !  READ  STRING  FROM  TAPE  I  ' 

READ  «i; I,D*C I > 

D*  <  I  >  =  D*  (I  >  C  1 , 3 1  9  3  S:C  HR*  ■:  1  3  > 

PRINTER  IS  16 

FOR  J=1  TO  200  !  BREAK  STRING  INTO  VARIABLE 

IF  J: 1  THEN  290 
F  C  1  ,'=POS  C  D*C  I  >  ,  CHR*  <1  3  >  ) 

L f  =  P 0 S  C  D#  C  I  >  C  1 ,  P  C  1)3,  CHR*  C  1  0>  ) 

L  =  LEN' D*C I ) C 1 , PC 1 )  1  > 

IF  Lf<>0  THEN  240 

IF  L<2  THEN  280 

0(1  )*VALCD#CI  )tl  ,  PCI  )  ]  > 

D  =  1 

GOTO  290 

IF  L< 3  THEN  280 

oci>=valcd*ci>c2,pci.<3> 

D=1 

GOTO  290 
Oc  1  >=-l 
Jl= J+l 

P  <  J 1  )  =ROS  <  D*  C  I )  [  P  C  J  > ♦ 1 3 , CHR# C 1 3 >  ) 

IF  P ( J 1 >  =  0  THEN  480 
PC J1 >=PcJl >+PC J) 

L=LENCD*C I >CP< J>  +  1 , PC  J1  >  3  > 

IF  L < 4  THEN  460 

Lf  =  POS  C  D*  CD,  CHR#  C  10  >  > 

IF  Lf <  >0  THEN  390 
G=  1 

GOTO  430 
G  =  2 

D* C I  >  C P  C  J  >+G, PC  J 1  >  3  =  TRI M* C D* c  I  )  t  F' C  J  >  +G ,  P C  J  1)3) 
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APPENDIX  A 

1 

4  id 

L=LEN< D#<  I  >  C  P  <  J  >  +  G, PC  J l >  3  > 

420 

IF  L  6  THEN  460 

4  3d 

0  C  J 1 >  »VRL  C  D# C I >  C  P C  J  >  +G, P  C  J 1 >  3  > 

44d 

Ii  =  H  -f  1 

450 

GOTO  470 

* 

460 

0  <  J 1  '=-1 

4  70 

GOTO  5O0 

480 

S*D 

, 

490 

GOTO  510 

1 

5dd 

NEXT  J  • 

NEXT  J 

| 

51  d 

8i  =  0 

52  d 

Ci*St + 1 -Ss 

53d 

C  f  =  3 1 +8 

1 

54d 

n=0 

i 

55d 

IF  C O < 1 >  =  - 1 >  AND  <S=0>  THEN  670 

1 

i 

56d 

F 0 R  K  —  C s  TO  C t  ! 

PACK  VARIABLES  FOR  OUTPUT  • k  > 

f 

5  7d 

L  =k-St  -t-Ss 

5Sd 

IF  0CL>=-1  THEN  600 

1 

590 

GOTO  640 

600 

L  =  L  +  1 

6 1  d 

Cf  =  C  f -1 

62d 

S  s  =  S  s  ♦  1 

6  3  O 

GOTO  580 

640 

F  K  :■  =  0  <  L  > 

6  5  d 

NEXT  K  ! 

NEXT  K 

66d 

St =  St +  S 

670 

NEXT  I  ' 

NEXT  1 

68d 

FOR  3  TO  32  ! 

OUTPUT  DATA  (ID 

6 '9d 

FOR  N=M  TO  1023+M  STEP  32 

700 

IMAGE  #,4D 

710 

PRINT  USING  700J FCN> 

7  20 

NEXT  N 

730 

PRINT  L I N< 1 > 

740 

NEXT  M  ! 

NEXT  M 

75d 

ON  END  #2  GOTO  790 

760 

CREATE  W2*,l,,9000 

77  0 

ASSIGN  #2  TO  W2# 

780 

PRINT  #2;FC*> 

790 

800 

|?hJT  W  "READ"  (Cont.) 

10  FOR  1=10  TO  50 

!  "COPY" 

20  I  #=VAL#  C  I  ) 

30  T  $  = " T 1  ID" 

!  FRAME  NUMBER 

> 

40  Plt=" : T15" 

!  TAPE  DRIVE 

50  P2$=" :F" 

60  F1#=T#&I#&P1* 

70  F2#=T#&I#S<P2# 

80  PRINT  F1#,F2# 

90  COPY  FI#  TO  F2# 

100  NEXT  I 

110  END 

!  DISC 

"COPY" 

1 

24 

i 

1 

i 

. 

_ y 
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10 

PRINT  "FRAME","  " , "  ",  "POSITION" 

!  "MRXAMP"  ! 

20 

FOR  1=10  TO  50 

30 

DIM  F<1 024) 

r 

40 

I*=VAL*< I ) 

50 

V#="T  1  :  F" 

60 

ASSIGN  #1  TO  V# 

70 

ON  END  #1  GOTO  100 

80 

READ  il;F<*> 

90 

Ma>:=0 

i 

100 

FOR  M= 1  TO  32 

!  OUTPUT  DATA  <M>  i 

110 

FOR  N=M  TO  1023+M  STEP  32 

1 

120 

IF  F < N ) >Max  THEN  140 

130 

GOTO  160 

140 

Max  =  F  <N) 

150 

Nax=N 

j 

1  60 

NEXT  N 

1 

170 

NEXT  M 

NEXT  M 

1  80 

PRINT  VS, "MAXIMUM  AMPL I TUDE " , Max , Nax 

i 

190 

PRINT  LINO) 

I 

200 

NEXT  I 

i 

210 

END 

1 

"MAXAMP" 


10  PRINTER  IS  16  !  "FREDIS" 

20  PRINT  PAGE;"  This  program  allows  the  user  to  enter  a  frame  of  data  and  ha 
"e  a" 

30  PRINT  "  frequency  distribution  computed  and  drawn." 

40  PRINT  LINC1); "  After  the  distribution,  a  beep  will  sound.  The  chart  will 


5£i  PRINT  "remain  on  the  screen  and  the  program  will  halt  until  you  press  CONT 

M 

60  PRINT  "This  will  give  you  time  to  digest  the  information  on  the  screen." 

70  PRINT  L I N < 1 > ; "  If  you  want  the  output  of  this  program  dumped  from  the  CRT 

II 


SO  PRINT  "to  the  internal  thermal  printer,  then  answer  Y  to  the  question  below 


90  PRINT  "If  you  don't  want  hardcopy  answer  N.  Then  press  CONT." 
100  INPUT  "Do  you  want  hardcopy  <YxN>?" , ft# 

110  IF  UPC#<A*>="Y"  THEN  Yhard 

130  IF  UPC# <  AS )  =  " N "  THEN  Nhard 

13©  BEEP 

140  GOTO  100 

150  Yhard:  Hardcopy=l 

160  GOTO  ISO 

170  Nhard:  Hardcopy=0 

130  N=224 

190  CALL  Dr i ver<N, Hardcopy > 

200  END 

210  SUB  Dr i ver < N, Hardc opy > 

220  RANDOMIZE 

2  30  DIM  fl(N),L$(N),Uniti$[30],VU024) 

"FREDIS" 
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I' 40  INPUT  "Which  frame  01  to  50?",  Fr* 

250  V*  =  "T1 1  V"&Fr*CV  :  F" 

2*0  ASSIGN  #1  TO  V* 

27©  READ  # 1 ; V < ♦ > 

280  Ma=  V  <  1  ) 

200  FOR  J=1  TO  1024 
300  Max=MAX ( Max , V ( J  > > 

310  NEXT  J 

320  MAT  A=ZER 

330  H i gh=0 

340  Lh i 9  h  =  4  O  9  6 

350  FOR  K= 1  TO  1024 

3*0  IF  (  V  (  K  >  >H  i  gh  >  AND  <  V  <  k  >  L h  i  9 h  >  THEN  380 
370  GOTO  390 
3 SO  H  1  g h  =  V (.  K  > 

390  NEXT  K 

40O  FOR  J 1 = 1  TO  1024 

418  IF  V ( J 1 >  *H 1 gh  THEN  43© 

420  GOTO  44© 

430  A  <  H  i gh ) =A  < H 1 gh  >  +  1 

44©  NEXT J  1 

450  IF  High*®  THEN  490 

4*0  Lhigh=High 

47©  High*© 

480  GOTO  350 
490  PRINTER  IS  0 
5©0  PRINT 
510  PRINT  V* 

52©  PRINT  "MAX  AMP" 

530  PRINT  Max 
54©  PRINTER  IS  16 

55©  FOR  1=20  TO  Max  STEP  2©  !  ENTER  THE  TITLES 

560  L$< I > =  VAL*  < I > 

570  NEXT  I 

5 8 ©  U n its**" N u i» b e r  o f  Pixels" 

590  T  i  1 1  e*=  "  Amp  li  t  tides " 

6©0  H  =  Max 

610  CALL  Ear <A<*> , L*<»> , N, Uni ts*,  Ti  1 1  e*> 

620  BEEP 
6  3  0  P  A IJ  S  E 

64©  IF  Hardcopy* 1  THEN  DUMP  GRAPHICS 
650  EXIT  GRAPHICS 
6  6  ©  S  U  E  E  X I T 

67©  SUB  Bar  <  A  <  *  )  ,  L*  <  *  ,  N ,  Un  i  t  s* ,  T  i  t  1  e* ) 

680  PLOTTER  IS  "GRAPHICS" 

690  GCLEAP 
70©  GRAPHICS 
71©  PEN  1 
720  Si.ini=0 


"FREDIS"  (Cont.) 
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7  50  Max=A'l> 

740  FOR  1=1  TO  H 
750  Surr»=Suift+fl< I > 

76U  Max=MRX(Max,  A<  I  >  > 

770  NEXT  I 
780  Ly=LGT<Max> 

790  SETGU 

300  LOCATE  15,120,40,95  !  Set  screen  boundaries  and  scale 

810  CLIP  0,123,40,100 
820  MOVE  61.5,97 
830  LORG  5 

840  LABEL  USING  850;Tit1e* 

850  IMAGE  #,K 

860  SCALE  0 ,  N+ .  5 , 0 ,  Max+Max  -*' 1 0 

870  Xtic=0  !  No  tic  marks  on  horizontal  axis 

880  Test  y  t i c  =FRACT  (.  Ly  >  +  < Ly<  0 ) 

390  Yt i c  =  10  - ( INTCLy >- 1 >*< 1  +  1 . 5*<  <Testyt i c >. 39794)  AND  <  Test yt i c  < . 69897  >  > +4* < < T e 
styt i c  >*.69897)  AND  <  Test  yt i c  <  =  . 87506  >>  +6 . 5*  <  Test  yt i c  >. 87506  >  > 

900  !  Compute  the  vertical  tic  marks 

910  CALL  LaxesCXt i c , Yt i c , 0, 0, 1 , 1 , 2, 0, N+. 5, 0, Max+Max/10> 

920  LORG  4 
930  FIXED  1 
940  FOR  1=0  TO  N 
950  MOVE  I - . 5, 0 

960  I  PLOT  0,A<I),-1  !  Draw  the  bar  for  the  ith  data  item 

970  I PLOT  1,0 
980  I  PLOT  0 , - A  < I > 

990  MOVE  I , A< I ) 

1000  NEXT  I 
1010  DEG 
1020  IMAGE  K 

1030  LDIR  90  !  Rotate  the  label  90  degrees 

1040  SETGU 
1050  LORG  5 
1060  MOVE  5,70 

1070  LABEL  USING  1020; Uni tsl  !  Write  the  display  units  being 

1080  !  used  along  the  vertical  axis 

1090  CLIP  0, 123,0, 100 

1100  SETUU 

1110  FUR  1=1  TO  N 


"FREDIS"  (Cont.) 
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1128  MOVE  I, -Max  20 
1136  LORG  8 
1 140  PEN  1 

1150  LDIR  90  !  Rotate  the  labels  90  degree 

1180  LABEL  USING  1020; L $(I> 

1170  NEXT  I 
1180  BEEP 
1190  SUBEXIT 

1200  SUB  Laxes<Xt-  i  c  ,  Yt  i  c  ,  Xorg,  Yorg,  Xmaj  ,  Ymaj  ,  M  i  rit  i  c  s  i  ze,  Xmi  n,  Xmax ,  Ym  i  n ,  Ymax 
1210  IF  <.  Xm  i  ri  >=Xmax  >  OR  ( Ymi  n>=Ymax>  THEN  SUBEXIT 
1 220  Xf udge  = . 82* < Xmax-Xm  i  n  > 

1 230  Yf  udge  = .  02*  <■.  Ymax-Y  min) 

1 240  AXES  Xt i c , Yt i c , Xorg , Yorg , Xmaj , Ymaj , M i m i c s i ze 
1250  DEG 

1260  IF  NOT  Xt i c  THEN  Labely 

1270  LabelxlIF  SGN<Xtic)*-l  THEN  Parx. 

1280  LDIR  90 
1290  LORG  8 
1300  GOTO  1330 
1310  Parx:  LDIR  0 
1320  LORG  6 

1330  FOR  I=Xorg  TO  Xmax  STEP  ABS<Xtic> 

1340  MOVE  I , Yor g-Yf udge 
1350  LABEL  USING  1360; I 
1360  IMAGE  #,K 
1370  NEXT  I 

1  380  Lately:  IF  NOT  Yt  i  c  THEN  SUBEXIT 

1390  IF  SGN<Yt i c  >■“!  THEN  Pary 

1400  LDIR  0 

1410  LORG  8 

1420  GOTO  1450 

1430  Pary!  LDIR  -90 

1440  LORG  6 

1450  FOR  I = Yorg  TO  Ymax  STEP  ABS< Yt i c  > 

1460  MOVE  Xorg-Xfudge,  I 
1470  LABEL  USING  1360;  I 
1480  NEXT  I 
1490  SUBEXIT 


"FREDIS"  (Cont.) 
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10 

DIM  F Cl 024) 

!  "PEAKS" 

26 

Ft=  "  T 1 1 V22 : F  " 

30 

L=  1 9 

40 

PRINTER  IS  0 

50 

PRINT  Ft,  "PEAKS  ABOVE", L 

60 

PRINT 

70 

PRINT  "ROW", "COLUMN", "POSITION". 

"AMPLITUDE" 

80 

PRINTER  IS  16 

80 

ASSIGN  #1  TO  Ft 

100 

ON  END  #1  GOTO  120 

1  18 

READ  #i;F<*> 

120 

FOR  M= 1  TO  32 

!  OUTPUT  DATA 

1  30 

C  =  0 

140 

FOR  N=M  TO  1023+M  STEP  32 

150 

C=C  + 1 

160 

IF  F  <  N  >  >  L  THEN  200 

170 

IMAGE  #, 4D, "  " 

180 

PRINT  USING  1 70; F<N> 

190 

GOTO  250 

280 

IMAGE  # , 4 D , “ * " 

210 

PRINT  USING  200; F(N) 

220 

PRINTER  IS  0 

230 

PRINT  M,C,N,F<N> 

240 

PRINTER  IS  16 

250 

NEXT  N 

260 

PRINT  L I N ( 1  ) 

270 

NEXT  M 

!  NEXT  M 

280 

PRINTER  IS  0 

290 

PRINT  LINO) 

300 

END 

"PEAKS” 
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10  DIM  F < 1 024 )  !  " DATAST " 

20  INPUT  "AMPLITUDES  HIGHER  THEN"  ,  L 
30  FOR  J=10  TO  50 

40  PRINT  LINO) 

50  PRINT  "FRAME  NUMBER ",J 

60  PRINT  LINO) 

70  S*=VAL*OJ> 

80  ASSIGN  #1  TO  "TllVlSn  " !  F" 

90  ON  END  #1  GOTO  110 

100  READ  #1JF<*) 

110  FOR  M=1  TO  32  !  OUTPUT  DATA  <M> 

120  FOR  N=M  TO  1023+M  STEP  32 

130  IF  FCN) >L  THEN  170 

140  IMAGE  # ,  4  D ,  "  " 

156  PRINT  USING  140; F<N> 

160  GOTO  190 

170  IMAGE  #, 4D, 

180  PRINT  USING  170JFCN) 

190  NEXT  N 

200  PRINT  LINO  ) 

210  NEXT  M  !  NEXT  M 

220  NEXT  J 

230  END 


"DATAST" 


30 
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Id  INPUT  "LEVEL" ,L 

2d  FOR  1=10  TO  50 

30  DIM  F  <  1 024  > 

40  I$*VAL*<I> 

50  «*"Tlir«$r:F" 

60  ASSIGN  #1  TO  V* 

70  ON  END  #1  GOTO  90 

80  READ  # 1 ; F  <  *  > 

90  GRAPHICS 

100  LINE  TYPE  1 

HO  SCALE  0,12,0,15 

120  FOR  M= 1 5  TO  25 

13©  C=0 

140  FOR  N=M  TO  1023+M  STEP  32 
150  C=C+1 

160  IF  F(NXL  THEN  200 
170  MOVE  M- 1 5 , C 

180  LABEL  F(N) 

190  IF  C=  1  3  THEN  210 

200  NEXT  N 

210  NEXT  M 

220  WAIT  3000 

230  GCLEAR 

240  SCALE  0, 13,0,230 

250  FOR  M= 19  TO  22  !  PLOT  COLUMNS 

260  WAIT  200 

270  C=0 

280  GCLEAR 

290  MOVE  0 , L  !  PLOT  LEVEL 

300  LINE  TYPE  4 

310  DRAW  13, L 


'  "FRAME" 


!  DATA  FROM  DISC 


!  LABEL  AMPLITUDE 


APPENDIX  B 


320  LINE  TYPE  1 
330  MOVE  1,150 
340  LABEL  "COLUMN" 

350  MOVE  0.0 

360  FOR  N=M  TO  1023+M  STEP  32 

3?0  C  =  C+ 1 

380  PLOT  C-l,F<N> 

390  PLOT  C ,  F  <  N  > 

400  NEXT  N 

410  WAIT  500 

420  MOVE  0,0 

430  NEXT  M 

440  CCLEAR  ,  PL0T  ROWS 

450  MOVE  0,0 

460  FOR  N 1=239  TO  335  STEP  32  ,  rL0T  LEVEL 

470  MOVE  0 ,  L 

480  LINE  TYPE  6 

4<?0  DRAW  13, L 

500  LINE  TYPE  1 

510  MOVE  1,150 

520  LABEL  "ROW" 

530  MOVE  0,0 
540  FOR  M»1  TO  13 
550  N=N 1 +M 
560  PLOT  M- 1 , F < H > 

570  PLOT  M,F<:N> 

580  NEXT  M 

590  MOVE  0,0 

600  WAIT  500 

610  GCLEAR 
620  NEXT  N 1 

630  NEXT  I 

640  END 


"FRAME"  (Cont.) 
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